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Thermal motion of the flux lines (FL) gives rise to fluctuating magnetic fields. These dy- 
■ namic fields couple to the nuclei in the sample and relax the nuclear spins. Based on a model of 
harmonic fluctuations, we provide a theoretical description of the nuclear spin-lattice relaxation 
(NSLR) process due to the fluctuating FLs in clean type-II superconductors. At low fields, the 
calculated longitudinal relaxation rate T{~ is enormously enhanced at temperatures just below 
^ \ T c . At intermediate fields, the resulting T{ 1 exhibits a peak structure as a function of tempera- 
' ture, which is eventually suppressed as the field is increased. The vibrational modes which have 
On components propagating along the FLs play an essential role in the T-f process. 
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Several recent experiments and theories have pointed to the importance of thermal fluctuations 
in the thermodynamics of type-II superconductors. Jl], || [| These fluctuations also have pronounced 
effects on the dynamic properties of the system. Recent measurements of the NSLR rate T{~ in 
the organic superconductor K-(ET) 2 Cu[N(CN) 2 ]Br (ET) by De Soto and Slichter et al. provide an 
excellent example of this. Q It was shown that there is an enormous enhancement of the relaxation 
rate of 1 H nuclei in the superconducting state relative to the normal state Korringa value. It is 
believed that the experimentally observed peak structure does not belong to the category of the 
Hebel-Slichter coherence peak, rather it is caused by the thermal motion of the fluxoid system. [|J 

In this letter, we attempt to provide a theoretical description of NSLR in clean type-II super- 
conductors induced by the fluctuating magnetic fields of the flux line lattice (FLL). In the presence 
of fluctuating FLs, the local magnetic field at nuclear sites varies in time, which may contribute to 
the NSLR. This mechanism becomes dominant in a system such as ET when relaxation through 
other channels (e.g. through the interaction with conduction electrons) is suppressed. We will 
confine our study to the region H <C H c2 , where the mean distance a between FLs is much 
larger than the radius of the vortex core £ (coherence length). Since we aim to give a general 
treatment of the phenomenon, we will, for simplicity, consider a model of continuous FLs based 
on the 3D London description. The case of layered superconductors with the external magnetic 
field perpendicular to the layers (chosen along the z-direction) could be handled similarly. 

We specify the axis of the j-th FL at time r by a 3D vector rj{r). The local field B(f, r) varies 
with time r only through the motion of the FL and is given by[[| 

b>» = & l/A) -ik-f^mZf^'™ ■ w 



where b(k) = B[zk\ — k±k z )/kj_(l + k 2 X 2 ), k = (k±,k z ) and k± = (k x ,k y ), A is the magnetic 
penetration depth, the line integral / dsj is along the j-th FL, V is the volume of the sample, 
and iV and L are the number and the length of FL respectively. To obtain the Fourier transform 
of B(f, t) we used V • B — 0. A similar expression for the layered system can be found in Refs. 
§1,0- 

It is seen clearly from Eq. (1) that the thermal motions of FLs lead to fluctuating magnetic 
fields. These fields are coupled to the nuclear moments in the sample. This coupling is described 
by the interaction TC(t) = —jn^J^n^ ' B(r n ,r), where 7 n is the nuclear gyromagnetic ratio, I is 
nuclear spin operator, and r n is the position of the n-th nucleus. We then have[^] 

^ = 27* E l(t Dl 2 / ^co8a;or(B a (f;0)B a (f;r)> , (2) 

il a=X>y J 27T 

where ujq = jnB, ( ) represents the thermal average of the FL configurations. Here we have taken 
an average over the relaxation rates of all the nuclei (denoted by the horizontal bar in Eq. (2)), 
which will be evaluated by spatial integration. Using Eq. (1) we recast Eq. (2) into 

1 = rtK^B) = 7 2 n j^- 3 J d^l^l^S&u*) , (3) 

where K(to , B) is the fluctuation spectral density function, S(k, too) = [S(k, cu ) + S(k, — uJq)], and 
S(k, ±u ) = J + °° ^e^_L. £ J dz . J ^(^J^W) (4) 

°° 33' 
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is the dynamic structure factor of the FLs, which is entirely determined by the FL structure and 
corresponding time evolution spectrum, without reference to any properties of the nuclei. We 
emphasize that our result (3) is quite general, applying even to the flux liquid and glass states. In 
the following we will focus our study on the FLL state. 

The problem of finding T^ 1 is now reduced to obtaining the dynamic structure factor S(k, &q) 
of the FLs, which is a fundamental dynamic quantity and of importance in many other applications 
(e.g. the inelastic scattering). Evaluation of S(k,u Q ) requires information about the dynamics of 
the FL fluctuations. Bulaevskii et al.@ have used the Langevin equations for the overdamped 
motion of pancakes in interpreting their T-f data in TLjE^CuOe, and found a monotonic decrease 
of Ti 1 with temperature. At this time the dynamics of the FLs are much less well understood than 
their static properties. A clean translational invariant superconductor is thus a valuable model 
system which can provide significant insight into the physics involved. In this case it is clear that 
the vortex motion is reversible. Once this fundamental case is understood, various complications 
(e.g. including damping of the vortex motion) can be introduced. 

For a clean superconductor, Fetter has derived a set of dynamic equations for the FLs by using 
the London model, which provides a basis for our study. |]] We assume that the deformation of 
a rectilinear FL is small so that the vortex axis can be specified by fi(r) = (Ri + Ui(z,T),z), 
where Ui(z,r) is the displacement of the i-th FL from its equilibrium position (R^ ,z). In the 
harmonic approximation to the total energy of the FLL, the motion of Ui(z,r) can be quantized 
by interpreting u xi (z) and u yi (z) as quantum mechanical conjugate variables.^] This leads to 

n xi (z,r) = ^E^-^^ + i/l^^* , (5) 



u yi (z,T) = ^E^[(« + H«- i V ! / iur + («-^)i,/1e*^ ) * ,(6) 

where ~R = h/2m, Q = Q + j^j^V + £ ; P is the mass density of the superelectrons, and the 
summation over q± is restricted to the first Brillouin zone. In the above equations, Q, a, rj, and 
£ are functions of (q±,q z ) and are defined in Eqs. (43) and (44) of Ref. ||. The corresponding 
Hamiltonian of the FLL reads H-fll = \Hq^q z h UJ (a q±qz a) q ^ qz + ct^ q±qz a q±qz ) , with the dispersion 
relation given by 

u 2 (q±, q z ) = (n + \^rf 2 v) 2 ~ (« 2 + f ) ■ (7) 
1 + q z z \ z 

Here cb\ ±qz and a^^ qz obey boson commutation relations, and represent the creation and annihilation 
operators for a single vibrational quantum (which will be referred to as a "fluxon" ) in the normal 
mode with wave vector (gl, q z ) and frequency u(q±, q z ). 

We will now derive S(k, u ) resulting from Fetter's dynamics. || As usual, we use the Bloch 
identity to write (^■u j (^) e -ik±-a j ,(z',r^ = e -2W(fc x ) e <fc^(*,o) ^.^V)^ where w(k±) = |([jfej 



Uj(z, 0)] ) is the Debye- Waller factor. The lowest order contribution arises from the one-fluxon 
process, in which a nucleus flips its spin and the FLL creates or annihilates a fluxon to conserve 
the total energy. Combining the Bloch identity and Eqs. (4)-(6), and keeping only the one-fluxon 
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term, we obtain 



S(k, uj ) 



2pn 



-2W(kj_) 



E E 7jP*£ + & k l ~ 2ak x k v }[n5(uj + w ) ^ 



<2_L<Jz if 



where K is the reciprocal lattice vector, Vt' = Q + Y+fp2^ ~ £> an< ^ n = ( a <Li3z a 9J-^ * s ^ e Bose- 



Einstein occupation factor of fluxons. 

Upon substituting Eq. (8) into Eq. (3), we obtain 



EE 



n 2 \ 2 



, coth 



_ , , c -2W«r t +K) 2k B T / g\ 



X 



[tt(q x + K x f + tt'(q y + K y f - 2a(q x + K x )(q y + K y )]8(u - u ) 



For arbitrary fields and temperatures the quantity in the summation is a complicated function of 
K and (q±,q z ), and the summation can only be evaluated numerically; however, in some special 
cases and with some restrictions, analytical results can be found. The simplest case is a single FL 
in a bulk sample (this also includes the case of well-separated FLs, i.e. a 3> A) . JTTJ] In this case 
the motion is entirely self- induced and only waves propagating along the FL (q± = 0) exist. Thus 
we have Q 7 = rj = £ = 0.||ln the long wave length limit the dispersion relation can be written as 
u)(q z ) ps ^r(hi y — 7 + f )(? 2 ,[@, 0, where 7 = 0.5772 is Euler's constant. We then obtain 
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The dashed curve 
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where Ei(x) is the exponential integral, |L3| q z0 is determined by to(q z0 ) 
in Fig. 1 shows T-f 1 as a function of t, where A = A / y/1 — t A , (u 2 )/X 2 

and t = T/T c . Here we have chosen uq = 5 MHz, 7„ = 2.675 x 10 A s~ l Gauss 
(3 = 0.1. 

As can be seen from Eq. (9), for arbitrary vortex density, T-f 1 is a sum of contributions from 
all normal modes. It is instructive to analyze the contributions arising from two special types of 
modes. First, we notice that the modes propagating perpendicular to the FL direction (q z = 0), 
which corresponds to the vibrations of the rigid FLs, do not contribute to T{~ , because there 
is no transverse field associated with these modes (as can be seen from Eq. (2), a fluctuating 
transverse field is required for longitudinal relaxation). This indicates the essential role played by 
the modes with q z 7^ in the Tf 1 process. Second, consider the vibrations propagating along the 
FLs (also known as the "helicon modes" || |12|). In the dilute limit these modes become identical 
to those of a single FL. At high fields (ao <C A), the dispersion relation of such modes reads 



w(q z ) PS ^(ln y ~ 7 + \ + 777 l n f )<? 2 i n the l° n g wave-length limit. In this case we also have 



K A ^> 1 for K 7^ 0. This allows us to keep the K = term only in Eq. (9) and to obtain the 



contribution ^- L - 



q_l_dq 1 



=0 from these so-called "helicon" modes IT 
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The NSLR rate T^ 1 arises, in general, from the contributions of both the "helicon" modes and 
those branches with q± ^ 0. But one should note that not all q± ^ branches contribute to the 
Ti process. This is seen clearly from Fig. 2, in which we have plotted uj(q±, q z )/too as a function 
of q z for a few values of q±. The branches with u(q±, q z = 0)/co>o > 1 do not contribute to T-f due 
to the 5-function in Eq. (9). For X H we have numerically evaluated Eq. (9) for ujq = 70 and 1000 
MHz (corresponding to B = 0.37 and 3.74 Tesla, respectively). The resulting T{~ 's are plotted 
as a function of temperature in Fig. 1 as the solid and dotted curves. Notably, a relaxation peak 
is found below T c . Note that the results shown in Fig. 1 only contain the contributions from the 
fluctuating FLL only. At very high fields, the relaxation peak could be absent if T C (B) becomes 
less than the temperature at the peak, or if the relaxation due to the FLL is so greatly suppressed 
that it is slower than that due to other sources, e.g. the interaction with the current carriers. 

We now address the behavior of T{~ at different fields in more detail. First of all, we mention 
that at typical NMR frequencies (from a few to several hundred MHz) the condition <C 1 
is satisfied for superconductors with T c ~ 10K. The low frequency region corresponds to fields 
of a few hundred Gauss. The FLs are well-separated, and are very soft in this case. T{~ can be 
described by Eq. (10) and from Fig. 1 we see that Tf 1 is enormously enhanced right below T C (B). 
Note that so far there seems to be no experimental data published in this regime. It would be 
interesting to carry out such measurements. When the field increases to clq ~ A, the interaction 
between the FLs begins to play a role. This interaction promotes the modes with q± ^ and 
stiffens the FLs. In the high frequency range (ao <C A), only the K = term gives important 
contributions to Tj -1 . The magnitude of Tf 1 is greatly suppressed in comparison with the single 
FL case, as can be seen in Fig. 1. A close examination of Eq. (9) suggests that this suppression 
comes from the decrease of the thermal factor coth 2 j^°, , the density of states of fluxons at uq, and 
the resonant wave number q z0 . These effects overcome the increase brought on by the prefactor 

in Eq. (9). One should note that the Fetter's vortex dynamics breaks down if the field is 
close to H c2 so that vortex fields overlap strongly. In this case the nonlocal effects (wave vector 
dependence of the elastic moduli) act to enhance the fluctuations, |2|, ||| and a Ginzburg-Landau 
type of approach becomes relevant. 

The fluctuation spectral density function K(u, B) contains detailed dynamic information of 
the system. It depends on both frequency uj and magnetic field B. In Fig. 3 we plot K(u, B) vs uj 
for a few B's of interest. At a fixed B, the function has a maximum at a frequency u m which shifts 
upward as the field increases. The field dependence of K(u, B) at a given frequency uj can also be 
seen from Fig. 3. At low frequencies, K(u, B) decreases with the field, which has been observed 
experimentally {|]. At high frequencies it increases with B, and at intermediate frequencies it has 
a peak structure as a function of B. These are very interesting predictions which can be tested 
experimentally by measuring the NSLR rates of different type of nuclei, or different transitions 
of a nucleus with I > 1/2.||15| The measurement of K(uj,B) provides a crucial test for a vortex 
dynamics. 

In summary, we have studied the NSLR caused by thermally fluctuating FLL in a clean type- 
II superconductor. Its applications for the vortex dynamics go beyond the NMR. [17] The NMR 
approach to the FLL state is quite diagnostic, and provides a sensitive probe to the dynamic 
fluctuations of the vortices. The simple harmonic fluctuations of the FLL seems to render a good 
description of the overall features of Tj obtained in ET. We found that the fluctuation effects on 
the T^ 1 process is most pronounced at low fields, where the FL is easy to "bend" and the thermal 
excitations of the fluxons with the frequency uj q is effective. In the above discussions, we have 
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ignored the effects of damping of the vortex motion. We expect that inclusion of weak damping 
would not change our results significantly, and that the theory should be directly applicable to 
systems with sufficiently large Hall angle[|16|. The harmonic and the overdamped vortex motions 
(this is believed to be the case in YBCO) represent two different important regimes. The latter 
also deserves further study. While at present there is still no well-accepted dynamic theory of the 
vortices, we believe that our calculations shed important light on this subject. Finally, we mention 
that the effects discussed in this paper should be more pronounced in layered superconductors since 
in that case the FLs are more flexible. 0] 
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grateful to CP. Slichter for communicating their experimental results prior to publication, and for 
drawing our attention to the calculation of spectral density function. Support by the ONR grant 
No. N00014-89-J-1157 is also gratefully acknowledged. 
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FIGURE CAPTIONS 



FIGURE 1. Ti 1 as a function of T/T c at c^o = 5 MHz (dashed curve, represents the case of a 
single FL), and at ojq = 70 (solid curve) and 1000 MHz (dotted curve). The unit of Tf 1 is chosen 
to be 87r 2p%2^ A 3 ; which is about 10~ 3 to 100 Hz for typical type-II superconductors. Note that we 

have plotted Tf 1 all the way up to T/T c = 1. In reality, the Tf 1 is only valid up to T C (B). Above 
T C (B), the system goes to the normal state. 

FIGURE 2. u(q z , q±)/co at cu = 1000 MHz and T/T c = 0.4 for q ± \ = 0, 0.15, 0.31, 0.46, 0.61, 0.76 
and 0.92 (from bottom to top). 

FIGURE 3. Fluctuation spectral density functions at B — 1000 (solid curve), 3000 (dotted curve), 
and 6000 Gauss (dashed curve). 
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